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[\]  The  multimodal  structure  of  the  baroclinic  tides  in  the 
Strait  of  Gibraltar  is  studied  using  observations  and 
numerical  simulations.  Observational  data  and  model 
results  revealed  the  presence  of  two  types  of  tidal  internal 
waves  generated  over  Camarinal  Sill  (CS),  One  propagates 
toward  the  Mediterranean  and  disintegrates  inlo  series  of 
nonlinear  short  internal  waves  with  amplitudes  of  50  m  and 
more.  The  second  type,  behind  the  first,  propagates  slower 
and  has  a  longer  wavelength.  The  vertical  structure  with 
both  upward  and  downward  displacements  of  isopycnals 
can  be  identified  as  a  manifestation  of  higher  baroclinic 
modes.  Analysis  of  the  empirical  orthogonal  functions  of 
the  ADCP  measurements  performed  over  CS  and  model 
lime  series  has  shown  that  the  second  baroclinic  mode 
predominates  in  the  second  type  of  internal  waves.  Its 
amplitude  can  reach  one  third  that  of  the  first  baroclinic 
mode  of  the  leading  waves  of  depression.  Citatum:  Vazquez, 
A..  N.  Stashchuk,  V,  Vlasenko,  M,  Bruno,  A.  lzquierdo,  and  P.  C. 
Gallacher  (2006),  Evidence  of  multimodal  structure  of  the 
baroclinic  tide  in  the  Strait  of  Gibraltar,  Geophys.  Res.  Lett,  33, 
L 1 76(15 ,  doi :  i  0. 1 029/2 006GL0268G6, 

1.  Introduction 

[:)  The  Slraii  of  Gibraltar  has  a  great  impact  on  the 
environment  of  the  Mediterranean  Sea  and  Atlantic  Ocean. 
It  can  make  an  important  contribution  to  tire  formation  of 
the  global  oceanic  circulation  by  supplying  nutrient-rich 
heavy  salty  waters  from  the  Mediterranean  Sea  to  the 
Atlantic  Ocean  via  the  Mediterranean  Out  How  (MO) 
[Rahmstorf,  1997;  Reid,  1979]. 

[}]  The  circulation  in  the  Strait  of  Gibraltar  is  roughly 
characterized  by  a  iwodayer  structure  with  Atlantic  waters 
(less  salty'  and  warmer)  flowing  into  the  Mediterranean  in 
the  upper  layer  and  the  Mediterranean  waters  (saltier  and 
colder)  flowing  out  below  [Armi  and  Farmer,  1988; 
Tsimplis,  2000].  The  density  of  the  final  MO  penetrating 
into  the  Atlantic  Ocean  and  forming  the  Mediterranean 
waters,  is  determined  by  the  mixing  and  enftainment  pro¬ 
cesses  at  small  scale  in  the  strait.  The  density  of  the  Atlantic 
water  inflowing  into  the  Alboran  Sea  from  the  strait,  must 
also  be  affected  by  these  processes. 

[a]  One  of  the  most  important  mechanisms  for  mixing  in 
the  Global  Ocean  is  internal  waves.  Particularly,  in  the  Strait 
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of  Gibraltar  high  amplitude  internal  waves  associated  with 
the  internal  tidal  dynamics  have  been  recognized  as  die 
main  contributors  to  the  mixing  between  Atlantic  and 
Mediterranean  layers  [Wesson  and  Gregg t  1988;  Bruno  et 
aL,  2002;  Macias  et  a/,T  2006].  This  strait  is  know  n  as  one 
of  the  most  famous  “hot  spots"  of  the  World  Ocean  where 
large  amplitude  internal  waves  are  regularly  observed  [see, 
e.g,T  Armi  and  Farmer,  1988;  Farmer  and  Armi t  1988; 
Wesson  and  Gregg.  1994;  Pettigrew  and  Hyde,  1990; 
Watson  and  Robinson,  1990;  Bray  et  aL,  1990;  Watson, 
1994;  Richez,  1994;  Pistek  and  La  Violet te,  1999;  Tsimplis, 
2000;  lzquierdo  et  at<%  2001;  Bruno  et  aL,  2002].  These 
waves  are  generated  by  tidal  flowr  interacting  with  Camari¬ 
nal  Sill  (CS)  which  is  located  near  the  western  entrance  to 
the  Strait  of  Gibraltar  (see  Figure  le),  According  to  the 
present  understanding,  an  interfacial  depression  is  formed 
over  the  western  flank  of  CS  when  the  tidal  curreni  is 
flowing  toward  the  Atlantic  Ocean  (outflow).  This  depres¬ 
sion,  being  arrested  over  CS  by  strong  tidal  flow,  is  released 
and  begins  to  move  eastward  when  the  outflowing  current 
slackens.  During  its  propagation  iowrards  the  Mediterranean 
Sea,  the  baroclinic  bore  disintegrates  into  a  series  of 
nonlinear  short  internal  waves.  As  the  intcrfacial  disturban¬ 
ces  arrive  to  Tarifa  Narrows  (TM)  (sec  Figure  le),  the  tidal 
current  starts  to  flow  eastward.  As  a  result  another  weaker 
baroclinic  bore  is  formed  over  the  eastern  flank  of  CS.  It 
propagates  westward  as  soon  as  tidal  tlux  weakens. 

[s]  Most  of  the  existing  works  devoted  to  the  study  of  the 
baroclinic  tidal  response  of  the  Strait  of  Gibraltar  have 
focused  on  the  analysis  of  these  waves,  which  in  fact 
demonstrate  an  overwhelming  first-mode  baroclinic  signal. 
However,  the  measurements  by  Armi  and  Farmer  [1988] 
and  Fanner  and  Armi  [1988],  as  well  as  observation  by 
Wesson  and  Gregg  [1988]  also  revealed  evidence  of  a 
higher  baroclinic  mode  signal,  which  has  not  been  described 
yet  by  any  theoretical  model.  In  this  paper  we  present  a  new 
experimental  data  set  which  show's  evidence  of  the  multi¬ 
modal  structure  of  baroclinic  tides  in  the  Strait  of  Gibraltar. 
These  data  are  compared  with  the  results  of  a  fully  nonlin¬ 
ear,  nonhydrostatic  numerical  model  Finally,  an  Empirical 
Orthogonal  Function  (EOF)  analysis  is  applied  to  the  model 
and  measured  velocity  data  in  order  to  quantify  the  modal 
structure  of  the  baroclinic  waves. 

2.  Experimental  Data 

[ft]  The  experimental  data  used  in  this  paper  were  col¬ 
lected  on  different  cruises  during  the  G1MIX  project.  Two 
typical  along-strait  temperature  sections  showing  the  multi¬ 
modal  structure  of  baroclinic  tides  over  CS  are  presented  in 
Figure  1.  They  were  built  on  the  basis  of  an  XBT  survey 
performed  from  the  vessel  moving  with  constant  velocity 
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Figure  1.  Temperature  transects  across  CS  obtained  on  (a)  1 1 .09.03  and  (c)  12. 1 1 .03,  fh  and  d)  the  associated  tidal  current 
predictions  over  CS,  and  (e)  bottom  topography  of  CS  and  the  map  of  the  Strait  of  Gibraltar  with  indication  of  the  transects 
for  Figure  [a  (red  line)  and  lb  (blue  line).  TN  denotes  die  Tariffa  Narrows. 


over  Camarinal  Sill.  Figure  la  was  constructed  from  34 
temperature  profiles  taken  between  12: IS  and  14:04  UTC 
on  the  1 1th  September  2003  with  a  mean  spatial  resolution 
of  270  m.  The  second  transect  (Figure  1c)  was  performed 
between  13:01  and  13:51  UTC  on  the  12th  November  2003. 
In  the  last  case  52  XBT  profiles  were  taken  with  a  mean 
spatial  resolution  of  170  m.  Figures  lb  and  Id  show  the 
prediction  of  the  along  strait  component  of  tidal  velocity 
over  CS  during  the  XBT  measurements.  This  prediction  was 
performed  on  the  basis  of  a  harmonic  analysis  using  38 
harmonic  components  of  velocity  at  45  m  depth  registered 
over  die  CS  from  October  of  1994  to  April  of  1995.  In  a 
near-bottom  layer  the  temporal  variability  of  die  currents  is 
also  characterized  by  tidal  periodicity  [Garcia  hafuente  et 
«/.,  2000]:  however,  ibe  absolute  value  of  velocity  here  is 
different  due  to  die  two  layer  structure  of  the  circulation  in 
the  strait  described  above, 

[7]  Although  the  temperature  transects  were  obtained 
from  different  cruises,  they  reveal  quite  similar  structure. 
Both  of  them  were  taken  approximately  I  hour  after  the 
maximum  outflow  of  the  current  when  internal  waves  have 
been  generated  over  CS  (find  thick  segments  in  Figures  lb 
and  Id),  The  areas  covered  by  Figures  la  and  1c  are  not 
identical  although  both  include  the  area  just  above  the  sill 


The  comparative  analysis  of  two  temperature  fields,  mea¬ 
sured  just  after  the  maximum  outllow,  allows  one  to 
conclude  that  qualitatively  they  are  quite  similar  even 
though  they  were  taken  two  month  apart.  At  the  cast  side 
of  CS  die  vertical  temperature  structure  is  almost  undis¬ 
turbed  {Figure  1c)  whereas  just  above  the  silFs  crest  the 
isotherms  arc  elevated  and  some  of  them  even  outcrop.  A 
complex  structure  of  the  temperature  field  was  located  to  the 
west  of  the  CS.  Superposition  of  several  barodinic  modes  at 
the  Icc  side  of  the  CS  is  evident  (find,  for  instance,  the 
fragment  A T  with  all  isotherms  displaced  downward,  and  the 
fragment  B,  with  isotherms  displaced  upward  and  down¬ 
ward,  these  can  be  identified  as  first  and  second-mode 
waves,  respectively).  This  pattern  in  general  coincides  with 
the  observations  obtained  in  November  2003  (Figure  Ic) 
where  similar  waves,  A  and  B,  can  be  identified.  It  can  be 
seen  that  the  vertical  scale  of  the  isotherms  displacements 
exceeds  100  m  behind  the  sill.  The  14°C  isotherm  is 
displaced  as  much  as  200  m  which,  to  the  best  of  our 
knowledge,  may  be  the  largest  internal  wave  amplitude  ever 
recorded. 

[a]  Although  we  are  using  temperature  as  a  tracer  to 
delineate  the  structure  of  the  internal  waves,  we  are  not 
implying  that  the  temperature  gradient  necessarily  generates 


2  of  6 


uim 


VAZQUEZ  ET  AL:  MULTIMODAL  STRUCTURE  IN  GIBRALTAR 


M7605 


the  restoring  force  for  the  waves.  Salinity  can  be  very 
important  in  determining  the  density  structure,  hence  the 
buoyancy,  in  the  strait.  However  the  isotherms  accurately 
reflect  the  displacement  of  the  water,  hence  the  amplitude  of 
the  internal  waves, 

[<>]  The  temperature  gradients  here  arc  weaker  in  com¬ 
parison  with  the  layer  of  seasonal  thermocline.  As  it  will  be 
shown  below  from  the  analysis  of  the  boundary  value 
problem,  the  maximum  wave  displacement  is  expected  at 
the  depth  of  about  100  m  which  roughly  coincides  with 
undisturbed  position  of  the  14°C  isotherm, 

[to]  Another  i mportan t  eonc l usion  from  t he  o bse r v at i o na I 
data  is  that  the  multi-modal  nature  of  baroclinic  tide  in  the 
area  of  CS  is  clearly  evident.  At  first  glance  these  higher 
mode  waves  (packet  B)  may  be  thought  to  be  a  mixing  area 
produced  by  the  first  baroclinic  mode  waves  [Wesson  and 
Gregg,  1988].  However,  this  structure  is  not  stationary,  as 
will  be  shown  below.  It  has  a  phase  speed,  which  is  different 
from  that  of  the  first  mode  wave,  it  propagates  eastward 
along  the  Strait  of  Gibraltar  revealing  that  it  is  a  progressive 
wave.  The  study  of  the  complex  multi-modal  nature  of  the 
phenomenon  is  the  topic  of  the  rest  of  this  paper. 

3.  Numerical  Modeling 

[i  1  ]  A  number  of  papers  have  been  published  during  last 
two  decades  on  numerical  modeling  of  the  baroclinic 
dynamics  of  the  Strait  of  Gibraltar,  All  attempts  to  model 
the  baroclinic  tidal  dynamics  in  the  strait,  from  the  first 
application  of  a  three-dimensional  ocean  circulation  model 
by  Hang  [1989]  until  now,  failed  to  describe  the  real 
characteristics  of  generated  waves  because  of  the  “hydro¬ 
static'**  approximation  used  in  most  of  them  (see  also  the 
very  recent  paper  by  Satin  mo  et  ai  [2004]).  These  difficul¬ 
ties  were  partly  overcome  in  works  by  Brandi  et  ai  [1996] 
and  Morozov  et  ai  [2002],  where  nonhydrostatie  effects 
were  incorporated  in  the  models  allowing  for  description  of 
the  structure  and  dynamics  of  a  baroclinic  tidal  bore. 
However,  the  detailed  structure  of  higher  mode  baroclinic 
modes  was  not  studied  in  these  papers,  because  the  model 
by  Brandt  et  al.  [1996]  was  two-layered  and  numerical  runs 
by  Morozov  et  al.  [2002]  were  performed  on  a  coarse 
horizontal  gnd  (200  m  grid  step). 

[12]  In  this  paper  a  fully  nonlinear  non-hydrostatic  o- 
coordinate  numerical  model  was  used  to  reproduce  the 
baroclinic  tide  generated  over  the  CS,  The  model  is  2.5D 
because  the  equations  are  laterally  averaged  to  account  for 
variations  in  the  cross-section  of  the  strait.  Details  of  the 
model  are  described  by  F7 asenko  et  ai  [2005],  The  calcu¬ 
lation  domain  extended  longitudinally  from  5.09°W  to 
6. 1 0° WT  covering  a  total  of  90  km.  The  horizontal  grid  step 
was  50  m  and  140  o-lcvels  were  used  in  the  vertical 
dimension.  The  bottom  topography  of  the  Strait  of  Gibral¬ 
tar,  shown  in  Figure  le,  was  used  in  the  model.  The 
topography  was  slightly  smoothed  to  prevent  small  scale 
noise  from  developing  in  die  simulations  (compare  the 
bottom  plots  in  Figures  I  and  2).  Initial  temperature  and 
salinity  profiles  were  taken  from  real  CTD  measurements. 
The  fluid  was  initially  at  rest 

[n]  Tli e  model  was  forced  by  the  M2  tidal  current  during 
three  semi-diurnal  tidal  cycles.  Figure  2  shows  the  model 
predicted  evolution  of  the  density  field  in  the  area  of  CS 


during  one  tidal  cycle.  The  run  starts  at  1  =  0  h,  when  the 
current  is  zero  at  CS  changing  from  inflow  (eastward)  to 
outflow'  (westward),  During  the  first  two  hours  of  the  tidal 
cycle  the  accelerating  barotropie  flow  produces  a  weak 
elevation  of  isopycnals  upstream  of  the  crest  of  CS  and  a 
deep  isopycnal  depression  (^50  m)  at  the  lee  side  of  it.  This 
wave  of  depression  is  arrested  by  the  accelerating  barotTopic 
flux  behind  die  obstacle  where  it  extracts  energy  from  the 
background  flow-  and  grows  in  amplitude  due  to  the 
permanent  feedback  with  the  barotropie  flow, 

[14]  In  fact,  the  initial  wave  at  this  stage  can  be  treated  as 
a  lee  wave  propagating  upstream  (the  inclined  line  e-c 
represents  the  line  of  equal  phase).  Normally  the  arrested 
lee  waves  possess  characteristics  of  both  first  and  higher 
baroclinic  inodes  (all  isopycnal  displacements  are  down¬ 
ward  in  section  a-a  and  isopycnal  displacements  are  both 
upward  and  downward  in  section  b-b).  The  first  and 
higher  modes  possess  different  phase  velocities  thus  they 
effectively  separate  during  propagation.  This  process  of 
separation  is  clearly  seen  in  Figure  2.  By  f  =  3  h  the  fastest 
first-mode  baroclinic  bore  separates  from  the  rest  of  the 
wave  field.  As  a  result  of  die  separation,  two  structures  of 
isopycnal  displacements  can  be  distinguished  at  time  /  =  3- 
4  h,  resembling  the  structures  observed  in  the  XBT  transects 
(Figure  2,  shaded  areas).  The  first-mode  baroclinic  bore  (A), 
in  which  the  isopycnals  are  depressed  toward  the  bottom  in 
the  w'hole  water  column,  is  located  just  over  the  Camarinal 
Sill.  The  second  structure  (wave  B)  represents  the  2-nd 
baroclinic  mode  with  isopycnals  displaced  upward  above 
the  pycnoclme  and  downward  below  it.  The  scope  of 
vertical  isopycnal  displacements  (from  trough  to  crest)  is 
comparable  to  that  found  in  the  temperature  measurements. 
When  the  outflow  current  slackens  (/  =  4  h),  bore  A  starts  to 
move  eastward  disintegrating  into  a  packet  of  first-mode 
nonlinear  internal  waves  while  the  rest  of  the  wave  field,  B, 
is  still  retained  at  the  west  side  of  the  CS,  It  is  released  at  /  - 
6  h  when  the  current  reverses  and  starts  to  generate  a  new 
system  of  waves  propagating  westward  (packet  C). 

[15]  Figure  2  shows  the  first  tidal  cycle.  Thus,  initial 
transient  processes  may  influence  the  results  shown  in 
Figure  2.  We  deliberately  chose  the  beginning  of  the  motion 
to  show  the  generation  mechanism  more  clearly  (from  the 
initial  cycle,  when  the  wave  field  is  not  contaminated  by 
other  processes).  Note,  however,  that  for  subsequent  tidal 
cycles  the  mechanism  of  the  weaves  generation  and  evolu¬ 
tion  qualitatively  remains  the  same,  although  the  structure 
of  the  waves  field  can  be  slightly  different  in  details. 

[ig]  The  mean  phase  speed  of  the  baroclinic  bore  (A)  in 
the  numerical  experiment  was  about  1.40  m  s  which  is  in 
agreement  with  the  values  computed  by  izquierdo  et  ai 
[2001]  and  Sannino  et  ai  [2004].  At  the  same  time  the 
phase  speed  of  the  structure  B  in  the  numerical  model  was 
about  0.60  m  s_l,  which  is  less  than  half  of  the  phase  speed 
of  the  wave  A.  Similar  values  of  the  propagation  speed  of 
second  baroclinic  mode  for  this  area  were  observed  by  Armi 
and  Farmer  [1988]. 

4,  EOF  Analysis 

[17]  An  empirical  orthogonal  function  (EOF)  analysis  of 
die  time  series  of  vertical  velocity  obtained  from  observa¬ 
tions  and  numerical  simulations  was  carried  out  to  gel  more 
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Figure  2.  Model  predicted  evolution  of  the  density  field  during  one  tidal  cycle.  The  arrows  show  the  relative  intensity  and 
direction  of  the  tidal  current.  Letters  A,  B  and  C  denote  particular  fragments  of  the  wave  field.  The  shaded  areas  at  the 
bottom  represent  the  smoothed  topography  of  CS. 


detailed  information  about  the  structure  of  the  generated 
waves.  The  vertical  structure  of  currents  over  CS  was 
measured  by  a  190  KHz  Acoustic  Doppler  Profiler  (ADP) 
during  25  hours  on  1 0th  November  2003,  The  temporal 


resolution  of  these  measurements  was  60  s.  and  the  bin  size 
was  equal  to  2,5  itl 

[i*]  The  structure  of  the  first  two  empirical  functions  of 
model  (Figure  3a)  and  ADP  (Figure  3c)  velocities  are  in 


Normalized  amplitude 
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Figure  3.  Characteristics  of  the  first  (gray  lines)  and  second  (black  lines)  modes  over  the  top  of  Camarmal  Sill 
Normalized  amplitudes  of  EOFs  calculated  from  model  (a)  and  ADP  (c)  vertical  velocities.  Appropriate  time  series  of 
vertical  velocity  amplitude  (at  the  depths  where  ihe  EOFs  have  maximum)  calculated  from  (b)  the  model  results  and  (d)  the 
experimental  data. 
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qualitative  agreement  and  resemble  the  first  and  second 
dynamical,  modes,  g,{r)  obtained  from  the  boundary  value 
problem 

gzz  +  (ATzJ/c)2#  =  0, 

*(0)  =*(-*/)  =  0. 

Here  Ar(z)  is  the  buoyancy  frequency,  c  is  the  phase  speed, 
and  H  is  the  water  depth.  The  amplitude  of  the  first  mode 
does  not  change  sign  over  the  water  column  and  therefore, 
all  the  isopyenals  oscillate  in  the  same  direction  with 
maximum  displacement  at  a  depth  of  about  100  m.  The 
second  EOF  changes  sign  at  a  depth  of  about  130  m 
indicating  that  isopyenals  are  moving  in  opposite  directions 
above  and  below  the  pycnocline.  This  behavior  is  in 
agreement  with  the  second  dynamical  mode. 

[i<>]  The  temporal  variability  of  the  maximum  vertical 
velocity  of  the  two  first  EOFs  over  CS  are  shown  in 
Figure  3d  for  the  observational  time  series  and  in  3b  for 
the  model  time  scries.  Time  t  =  0  corresponds  to  the 
beginning  of  a  tidal  cycle  when  the  background  tidal  flow 
is  zero.  The  second  tidal  cycle  (from  /  =  T)  was  used  for  the 
EOF  analysis  of  the  numerical  data  to  insure  that  the 
"spin-up”  of  the  model  was  completed  and  no  transient 
processes  are  present  in  Figure  3  b.  It  should  be  noted  that 
EOF  analysis  was  carried  out  after  the  elimination  of  the 
temporal  mean  current  from  the  vertical  velocity  scries, 

[20]  Both  numerical  and  experimental  curves  reveal  ev¬ 
idence  of  quite  strong  first-  and  second-mode  signals.  An 
analysis  of  these  time  series  shows  that  the  maximum  signal 
in  both  cases  takes  place  3-6  h  after  the  beginning  of  the 
tidal  cycle.  As  was  demonstrated  above  by  the  numerical 
modeling,  the  generation  of  a  strong  depression  of  the 
pycnocline  above  CS,  its  separation  into  a  first  mode 
baroclinic  bore  (packet  A)  and  into  the  multi-modal  wave 
structure  (packet  B)  takes  place  during  this  time  span, 

[21]  Some  discrepancies  between  simulated  and  observed 
data  can  be  explained  in  terms  of  the  assumptions  made  in  the 
model.  For  instance,  they  could  be  caused  by  the  manifesta¬ 
tion  of  3D  effects,  which  were  not  included  in  the  model 
(except  for  the  variable  cross-seciion),  or  they  could  be 
caused  by  disregarding  the  mean  llow  through  the  strait. 
Another  possibility  could  be  the  use  of  smoothed  bathymetry 
in  the  model.  However,  apart  from  these  discrepancies,  one 
can  conclude  that  EOF  analysis  of  both  theoretical  and 
experimental  scries  revealed  qualitatively  similar  general 
patterns  with  evidence  of  the  presence  of  strong  second- 
mode  signal.  The  vertical  structure  of  the  first  two  EOFs  is 
very  similar  to  those  obtained  from  the  boundary  value 
problem.  In  conjunction  with  the  data  presented  above  from 
the  XBT  surveys  and  from  the  model  results,  the  conclusions 
derived  from  the  EOF  analysis  arc  strong  evidence  of  the 
multimodal  nature  of  baroclinic  tide  in  the  Strait  of  Gibraltar. 

5,  Discussion  and  Summary 

[22]  A  strong  semi-diurnal  tidal  forcing  plays  a  funda¬ 
mental  role  in  the  dynamics  of  the  Strait  of  Gibraltar.  This 
oscillatory  flow  interacts  with  Camarinal  Sill  creating  inter¬ 
nal  waves  which  are  the  main  source  of  energy  for  mixing 
and  entrainment  processes  between  the  Atlantic  inflow  and 


the  MO,  The  MO  spreads  into  the  eastern  North-Allantic 
forming  a  warm  and  salty  tongue  at  intermediate  depths, 
which  is  thought  to  contribute  to  the  preconditioning  of  North 
Atlantic  Deep  Water  formation,  and  hence  to  the  global 
thennohaline  circulation  [Reid,  1979],  In  addition,  the  mixing 
processes  induced  by  the  internal  waves  generated  at  the  CS 
of  the  Strait  of  Gibraltar,  modify  the  density  of  the  Atlantic 
water  jet  entering  the  Alboran  Sea  [Macias  et  al 2006]. 

[23]  The  analysis  of  observational  data  collected  in  the 
area  of  CS  [Bnmo  et  al. ,  2002]  and  the  results  of  numerical 
modeling  both  revealed  two  types  of  l urge -amplitude  inter¬ 
nal  waves  with  comparable  amplitudes  generated  by  tidal 
cuircnt.  The  first  one,  a  packet  of  large  amplitude  internal 
solitary  waves,  is  the  result  of  disintegration  of  the  well 
documented  first-mode  baroclinic  bore  generated  by  tidal 
flow  over  CS  and  propagating  towards  the  Mediterranean. 
The  second  system  of  waves  (behind  the  first)  propagates 
slower  and  has  relatively  larger  length  scale.  This  system 
reveals  the  characteristics  of  quasi-steady  higher  baroclinic 
mode  internal  waves  trapped  by  the  flood  flux  just  over  the 
sill  during  the  flood  phase,  lasting  for  4  h  (fragment  B  in 
Figure  2). 

[24]  From  the  observations  this  broadening  or  pycnocline 
might  appear  to  be  a  mixing  area  produced  by  the  first 
baroclinic  mode  waves  in  the  areas  of  strong  shear.  How¬ 
ever,  the  numerical  modeling  has  shown  that  this  structure  is 
not  stationary  over  the  whole  tidal  cycle;  it  has  a  phase 
speed,  and  propagates  eastward  along  the  Strait  of  Gibraltar 
resembling  a  progressive  wave  behavior.  The  vertical  slme- 
ture  of  isopycnal  displacement  (Figure  3)  along  with  the 
vertical  structure  of  horizontal  velocity  obtained  from  the 
numerical  runs  (not  shown  here)  clearly  indicates  lhat  this 
disturbance  can  definitely  be  treated  as  progressive  second 
mode  wave.  Thus,  the  experimental  data  presented  here 
along  wilh  their  theoretical  interpretation  can  be  considered 
as  a  strong  evidence  of  the  complex  multimodal  structure  of 
baroclinic  tide  in  the  Strait  of  Gibraltar.  Not  only  the  well- 
known  first- in  ode  internal  waves  but  also  higher  mode 
baroclinic  disturbances  with  comparable  amplitude  should 
be  taken  into  account  in  a  general  balance  of  mixing 
processes  in  the  strait. 
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